D- Algebra Structure of Topological Insulators 
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In the quantum Hall effect, the density operators at different wave-vectors generally do not com- 
mute and give rise to the Girvin MacDonald Plazmann (GMP) algebra with important consequences 
such as ground-state center of mass degeneracy at fractional filling fraction, and Wi+oo symmetry 
of the filled Landau levels. We show that the natural generalization of the GMP algebra to higher 
dimensional topological insulators involves the concept of a D-algebra formed by using the fully 
anti-symmetric tensor in D-dimensions. For insulators in even dimensional space, the D-algebra is 
isotropic and closes for the case of constant non- Abelian F(k) A F(k) ... A F(k) connection (D-Berry 
curvature), and its structure factors are proportional to the D/2-Chern number. In odd dimensions, 
the algebra is not isotropic, contains the weak topological insulator index (layers of the topological 
insulator in one less dimension) and does not contain the Chern-Simons 9 form (FA A — 2/3AAAAA 
in 3 dimensions). The Chern-Simons form appears in a certain combination of the parallel transport 
and simple translation operator which is not an algebra. The possible relation to D-dimensional 
volume preserving diffeomorphisms and parallel transport of extended objects is also discussed. 

PACS numbers: 74.20. Mn, 74.20. Rp, 74.25. Jb, 74.72. Jb 



INTRODUCTION 



Fractional topological insulators are the strongly cor- 
related states that may appear when a narrow bandwidth 
bulk band of a topological insulator [l|, Q is fractionally 
filled and subject to strong interactions. The plethora of 
new experimental facts and theoretical ideas discovered 
in the non-interacting topological insulators suggests that 
their interacting fractional counterparts will also exhibit 
new physical properties of topological phases, especially 
in space dimensions higher than 2. Evidence for their ex- 
istence has been provided in a series of analytical and nu- 
merical works in two-dimensional Chern insulators [3]- 11 
and time-reversal invariant topological insulators |12l Il3 
although more investigation is needed to determine the 
nature of the ground-state in the time-reversal invariant 
systems 14 1 . 



The most studied fractional topological insulator, the 
Fractional Chern Insulator (FCI), exhibits many of the 
properties of its Fractional Quantum Hall counterpart if 
the Berry curvature is smooth enough in the Brillouin 
zone (BZ). As evidence for the FCI state, it was analyt- 
ically shown in [n| that the counting of the quasihole 
excitations of the FCI state (in both the energy and the 
entanglement spectrum) is identical to the counting of 
quasiholes in a FQH state plus a lattice folding from the 
many-body FQH Brillouin zone to that of the FCI. The 
existence of a many-body Brillouin zone and the center 
of mass degeneracy of q at filling p/q in the FQH, as 
well as the mapping of the counting of states from the 
FQH to the FCI, is a direct consequence of the Girvin- 
Plazmann-MacDonald (GMP) algebra[l6j exhibited by 
the FQH translational operators. This algebra has a sim- 
ple interpretation in that the commutator of two trans- 



lation operators is non-zero and roughly proportional to 
the magnetic flux thru the area subtended by the trans- 
lations. The Chern insulator in 2-dimensions satisfied 
a similar algebra for parallel transport translational op- 
erators in the Brillouin zone 17]. This algebra has far 
reaching consequences: it is identical to the algebra of 
area-preserving diffeomorphisms, thereby providing for 
an explanation of the edge modes of an integer quantum 
Hall liquid as shape deformations of the liquid droplet. 
It allows for the construction of nontrivial many-body 
symmetry operators of the Hilbert space, it provides for 
a center of mass degeneracy (exact in the FQH but ap- 
proximate in the FCI), and is related to the Hall viscosity, 
<7 4 form factor, as well as the edge dipole-moment (l8| . 

In view of these consequences, the generalization of 
the GMP algebra to higher dimensions is desirable. In 
this paper, we present such a generalization. Realizing 
that the commutator of the projected densities p q in the 
" valence" bands of a topological insulators for two wave- 
vectors is a simple contraction with the antisymmetric 
tensor in two dimensions, we generalize this structure in 
D spatial dimensions to a £>-commutator by contract- 
ing with the antisymmetric tensor in D-dimensions. If 
this commutator is closed, the relation is called a D- 
algebra. We find that for topological insulators in even 
dimensions, the commutator is closed, and the algebra 
is isotropic, under a condition similar to that of the ex- 



istence of the GMP in the 2-D Chern insulator [1 71]: the 
F(k) A F(k) ... A F(k) object (which we call the D-th 
Berry curvature since F(k) is the Berry curvature in 
the BZ) connection must be smooth and ideally propor- 
tional to the identity matrix. The structure factors of 
the algebra are then proportional to the D/2'th Chern 
number. In odd space dimensions, where we would ide- 
ally expect the Chern-Simons form (which in 3 D reads 
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F A A - 2/3^4 A A A A), we find that this term does not 
appear in the algebra, which is now anisotropic and con- 
tains layers of (D — l)-dimensional topological insulators. 
While part of the non-Abelian Chern-Simons form (and 
the entire Abelian Chern-Simons form) appears in the 
odd--D-commutator, once we try to close the commutator 
to obtain a closed algebra, the Chern-Simons term van- 
ishes. We can, however obtain the Chern-Simons term 
as a combination of parallel transport and simple trans- 
lation operators, but they do not form an algebra. We 
close with comments and conjectures about the volume 
preserving diffeomorphisms and parallel transport of ex- 
tended objects that the D-algebra structure suggests. 

FLAT-BAND HAMILTONIAN AND PROJECTED 
DENSITY OPERATORS 

We start by fixing notations and recalling some well 
known results about band structure and projected den- 
sity operators for topological insulators. We consider a N 
band topological insulator described by a translationally 
invariant Hamiltonian, and we work on a D dimensional 
lattice with periodic boundary conditions. The Hamilto- 
nian reads as 

H = Aeh-liaMfi CO 

where a, (3 = 1, ...,N contain orbital and spin indices. 
We work in units where the lattice spacing is unity. The 
system being explicitly translation invariant, it decouples 
in momentum space 

H = E c LV/3(k)c k/3 (2) 

k,a,P 

keeping in mind that momenta are restricted to the first 
Brillouin zone. The Bloch Hamiltonian matrix h a @ (k) = 
Y^ r e~ ik ' T h r> oc,fi can be diagonalized Y^,p p = 

E n (k)u^ a , where the eigenvectors u£ ^ are chosen to 
form an orthonormal basis. Consequently the Bloch 
Hamiltonian can be separated into normal mode oper- 
ators 7^ at momentum k of the band n 

H = J2En(^^ (3) 

k . n 

where the normal modes can be written as a matrix rota- 
tion of the original electron operators 7^ = u^*gCu,p- 
We then define the projector into the band n at momen- 
tum k by 

P, i ,k = 7 k it |0)(0| 7 ^ (4) 

From now on, we will consider the physics of the (possibly 
fractionally) occupied bands and look only at projectors 



into the occupied bands. The projection operator in the 
occupied bands is 

p = J2 p ™± ( 5 ) 

where the band index n = 1, . . . , N occ ranges over all 
occupied bands. The projected density operator in the 
occupied bands of a topological insulator is the following 
operator: 

<°q = P E ^S^P = £ WK +q )i k nt |0> (0| 7k% 

j k,n,m 

(6) 

where n, m ranges over the set of occupied bands. A 
simple calculation reveals that: 

[H, Pq \ = ]T (S„(k)-i? m (k+q))«J< +q )7 k it |0) (0| 7 kV 

k.n,m 

(7) 

It follows immediately that the projected density opera- 
tors commute with the Hamiltonian if E n (k) = E m (k+q) 
for all occupied bands and all momenta. This is called 
the flat-band Hamiltonian : 

Hfb = - J2 p »i,k+ P » 2 .k- (8) 

n!,E ni <fj, n 2 ,E n2 >fj, 

Fractional topological insulators are usually constructed 
and observed in models with fractionally filled bands 
whose bandwidth is very small, such that interactions 
and not the kinetic energy dominate the physics. The 
ideal example of such an insulator is precisely the flat 
band one-body deformation (jHJ). 

We stress that projected density operators are an exact 
symmetry of the flat-band Hamiltonian © , to which the 
true one body Hamiltonian (JTJ is adiabatically connected. 
As a matter of fact, this also holds for any projected op- 
erator that does not mix occupied and unoccupied bands 
(in particular any projected operator) . 

THE TWO DIMENSIONAL CASE 

Before moving to higher dimensions, we quickly review 
what is known about the algebra of projected density op- 
erators in two dimensions, with an emphasis to its main 
characteristics. We focus on the appearance of the Chern 
number in the algebra and on the link between projected 
densities and parallel transport in the background of the 
Berry curvature. 

Density algebra and first Chern number 

At long wavelength qi, q 2 -> 0, the reference finds 
the following commutation relation : 
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■ U V 



k,n,m 



^;r(k)7k f io) <o| 7 r+c 



(9) 

where the Einstein summation convention after repeating 
indices is assumed. F^ v — d^A„ — d u A^ — i[A^, A v ] is the 
non-Abelian Berry field strength in the Brillouin zone, 
while the vector potential is A^ m (k) = i{u^.\d k \u^). 
This result holds in any dimension. 

In two dimensions i^„(k) = £?(k)e^„, and its integral 
over the whole Brillouin zone yields the first Chern num- 
ber 



C x = ±- f d 2 fce^Tr(F^(k)) 
47r Jbz 



(10) 



The commutator of two densities has to be non-zero in 
a nontrivial Chern insulator. That is so because the 
Chern number G\ of the bidimensional insulator can be 
expressed as a trace over the Brillouin zone of the density 
commutator 



Tr ([Pqi, Pq 2 ]P- 



qi-q2 



))^o^-(qiAq 2 )C 1 (11) 



where qi A q 2 = 9 2 ■ I n the continuum limit of the 
Quantum Hall effect, the projected density algebra of 
the Lowest Landau Level is called the CMP algebra. Its 
generators are the generators of the area-preserving dif- 
fcomorphisms in two-dimensions. This result is recovered 
for bidimensional topological insulators with an Abelian 
U(l) uniform Berry curvature, in the long wavelength 



limit. As pointed out in [17j (see also [19j and [15j]), if 



the local Berry curvature can be replaced by its average 



F xy {k) = B 



Ci 
2tt 



then 



[Pqi , Pq 2 ] = ~iB A Pqi+q 2 



(12) 



(13) 



Note that qi Aq 2 is the area enclosed in the parallelogram 
delimited by qi and q 2 . This algebra is nothing but 
the bidimensional Aharonov-Bohm effect in momentum 
space, in the background of the "magnetic field" F xy — 
B. 

This Abelian treatment applies to two-band models 
(insulators with one band below and above the gap) or 
to many-band insulators where the non-Abelian compo- 
nents of the field strength can be neglected (up to a pref- 
actor N occ ). 

Expanding the projected densities at long wave- vectors 
as p q = 1 + iq ■ R + 0(q 2 ), the algebra of the guiding 
center is recovered : 



[Rn, Rv] = iBe^y — — G\ tnu 
Ztt 



(14) 



The Chern number quantifies the non-commutativity of 
the guiding center operators. 

We remark that in a two-band insulator, it is impossi- 
ble to have a constant Berry curvature due to the no-hair 
theorem although this seems possible in insulators 
with four or more bands I20i . 



Regularized density operators 

Since projected density operators commute with the 
fiat-band Hamiltonian (|SJ|, it would seem that they are 
the generators of a proper symmetry group of the sys- 
tem. However this is not quite true, as they suffer from 
a serious deficiency. Because of the projection, they are 
not unitary. The density operator translates states in 
momentum space but does not keep their norm: 



p q |n,k> = £«_ q | 



(15) 



where |n,k) = 7]? |0). It is possible to replace the pro- 
jected density operator p q by a unitary operator p q , 
while not spoiling the long wavelength behavior from 
Eq.@. For a uniform Abelian Berry curvature, the an- 
swer is quite straightforward, and is simply the expo- 
nentiation of the guiding center operator. Doing so, one 
recovers the GMP algebra 



[Pqi j Pq 2 



qi A q 2 ■ 



„. • / D qi Aq 2 \ 
2i sin yB J 



Pqi- 



(16) 



More generally, for a non-Abelian and non-uniform Berry 
field strength, the answer is parallel transport in the 
background of the Berry gauge potential A^(k) : 

Pq= £ he-^^') 7 k nt |0)(0| 7k "V q . (17) 

k;n,m 

In the Abelian case this result was pointed out in [l7j . 
Note that this operator also commutes with the Flat- 
Band Hamiltonian, and at small momenta we recover the 
projected density operator p q = p q + 0(q 2 ). 



DENSITY ALGEBRA IN EVEN-SPACE 
DIMENSIONS 

The density commutator is natural in two dimensions. 
In higher space dimension D > 2, the commutator alge- 
bra Eq.((n]) reveals whether a two dimensional quantum 
Hall effect exists on planes of the D-dimensional space. 
For example, in D = 3, the existence of a commutator 
algebra Eq. (jH) for qi , q 2 in the xy planes for each z po- 
sition (momentum) reveals the existence or absence of a 
2-dimensional quantum Hall effect on those planes, triv- 
ially translated in the z direction. It is apparent then that 
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in D > 2 dimension the commutator algebra Eq. (j9)) can- 
not be isotropic. In order to find an isotropic algebraic 
structure in higher dimensions, we must look somewhere 
else. We first realize that the commutator [p qi ,/9 q2 ] is, 
in two dimensions, simply a re-writing of the operators 
Ea^PqaPq^ • This object can now be generalized to any 
dimension. In D space dimensions, it is then suggestive 
to look at the operator: 



Pq ai ) Pla 2 i ' ' ' i P<ia D — e a 1 a 2 - a D Pq al Pq a2 ' - ' Pq a£) 



(18) 



where e ai ,a 2 ,...,a D is the totally antisymmetric tensor in 
D-dimensions. and a = 1 . . . D. These generalized com- 
mutators are called D-commutators. We will now com- 
pute this object in the longe wavelength limit and find it 
is closed, thereby generating a D-algebra. 



D-commutator and Berry curvature density 

The density algebra in even space dimensions is simpler 
to obtain than in odd-space dimensions for reasons that 
will become apparent. In even space dimensions we have 
the Chern-insulator (QH)-classes, so we anticipate that 
the algebra closes. We first re-express the D-commutator 
as a product of 2-commutators: 

[Pqn Pq2 ! ' ' ' ) Pqo] = ^ ^ £a 1 a-2...a D - 1 a D 

[Paa 1 ,Pi a2 ]---[p<i« D _ 1 ,P<i aD ] (19) 

Using the long wavelength two dimensional algebra 
Eq.® and working at order q° we obtained : 

[Pqi.Pqa.-" »Pq D ] = H/2) D/2 (qi A q 2 A • • • A q D ) X 
(F(k) A • • • A F(k)) nm 7 »t |0) (0| 7^ qi+ ... +qD 

k . n , m 

(20) 

where we used the identity: 

z aia2 -.. aD <f a \cf a l ■ ■ ■ < D D = (qi A q 2 A • • • A q D ) e ^~^ 

(21) 

This equation is the D-dimensional analogue of Eq.(|9]). 
In the D-commutator appear the matrix 

F(k) A • • • A F(k) = e^-^F^ 2 (k) ■ • -F Mfl _ lMD (k) 

(22) 

which is the D/2'th Berry curvature density of the 
D/2'th Chern number: 

Cn,. = {D/m » /2 ^ )D/2 J d D kTr (F(k) A - A F(k)) 

(23) 



For even dimensional topological insulators, the D/2'th 
Chern number can be expressed as the the trace over the 
D-commutator of the projected density operator: 

([PquPq2J - ' ' )Pqr>]P-(qi+...+q.D)) ~q^0 
(2 TTl) D / 2 ( ^ /2)! A ^ A ' ' ' A ^ ° D/2 ( ^ 

This is the exact analog of the bidimensional relation 

CD. 



GMP in higher dimensions and D-algebra 

It is possible to obtain an analog of the GMP algebra 
in D-dimensions. As for topological insulators in two 
dimensions, this algebra holds when the Berry density 
F(k) A • • • A F(k) is uniform in the Brillouin zone, and 
proportional to the identity matrix. Then the projected 
density operators algebra becomes a D-algebra, in the 
long wavelength limit: 



[Pqi ! Pq 2 i' " i Pq 



(D/2)l- 



1 



C 



D/2 



(2tti) d / 2 N occ 
(qi Aq 2 A ■ • • Aq D )p qi +...+q D 



(25) 



and we recover a D-algebra. It is very tempting to ex- 
pand the projected densities as p q = l+iq-R+0(q 2 ). As 
a consequence of Eq. ([25|) the "guiding center" operator 
R also form a D-algebra: 



2tt 



D/2 n 

^D/2 



Jv occ 



(26) 



This algebra is most easily understood in the continuum 
limit. From Eq. (ll5[) the guiding center operators are sim- 
ply the covariant derivative with the Berry potential in 
momentum space: 



fl„(k) = -t (e fcM -a M (k)) 



(27) 



Using the relation = iF^, it is straightforward 

to obtain their D-commutator 



[Ji Ml (k),i? M2 (k),.- - ,#^(k)] = 
(i/2) D / 2 F(k)A---AF(k)e Ml ... MD 



(28) 



This elementary derivation in the continuum is not 
plagued by the limitations of the projected density oper- 
ators algebra, as there is no need to suppose the Berry 
density F(k) A • • • A F(k) to be uniform or proportional 
to the identity. 

This D-algebra structure may be understood in two 
ways. On one hand, as was pointed out in [TBI ], the pro- 
jected position operators can be expressed in terms of 
the projected density operators. Therefore an immediate 
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interpretation of (|26f is the non commutativity of the co- 
ordinates of particles projected to the occupied bands of 
a topological insulator. This is the D-dimensional analog 
of Ea.([T3)l for the Quantum Hall Effect. 

On the other hand, the GMP algebra also describes 
a bidimensional Aharonov-Bohm effect: the projected 
density operators implement parallel transport of point 
like objects in the background of the Berry curvature 
F. In higher dimensions, an Aharonov-Bohm effect with 
respect to the Z?-form F A • • • A F requires parallel trans- 
port of higher dimensional objects. Indeed, the natural 
objects that can couple to a D-form are D—2 dimensional 



membranes [21|. Interestingly, the classical limit of the 



22], which 



D-commutator is the Nambu-Poisson bracket 
is a natural setup to describe the dynamics of classical 
membranes [23j |. The appearance of extended objects in 
the field theory description of topological insulators in di- 
mensions greater than three is also expected from the BF 
proposal of [3]. We conjecture that the algebra (|28| is 
related to an Aharonov-Bohm effect involving extended 
excitations (membranes) coupled to the Berry curvature 
F AF A - • • F. However, unlike in two dimensions, it is not 
clear how to interpret the projected density operators as 
an implementation of membrane parallel transport. 



DENSITY ALGEBRA IN ODD-SPACE 
DIMENSIONS 

Pursuing the same strategy in odd dimensions leads to 
an impasse. The topological invariant in odd dimensions 
is defined as the integral over the Brillouin zone of a 
Chern-Simons form. For instance in three dimensions 
the Z2 topological invariant is given by 



= — 



1 



2tt 32tt 2 



d 3 k e^Tr 



F A 



~A.n A< 
3 M 



A,. 



(29) 

Defined for all odd dimensions, a characteristic feature of 
Chern-Simons form is that their integral is not necessarily 
invariant under large gauge transformations. However 
the variation has to be an integer (25[. Contrary to the 
even dimensional Chern numbers, the odd dimensional 
Z2 topological invariant is only gauge invariant modulo 
integers. Trying to obtain P3 through the gauge invariant 
trace 



(30) 



T r ([Pqi 1 Pq2 1 Pqa]/ 3 — qi— q2— q3 I 

is doomed to fail. A simple relation like Eq. ([23|) is ruled 
out in odd dimensions. Moreover £>-commutators in odd 
dimensions are known .26] be more problematic than 
their even dimensional counterpart. For instance while 
even commutators involving the identity matrix do van- 
ish, this is no longer the case for odd commutators. This 
is most easily seen in 3 dimensions: 



[A,B,1] = [A,B]?0 



(31) 



Consequently, when expanding the 3-commutator of a 
projected density operator p q = 1 + iq • R + 0(q 2 ), the 
lowest order contribution is of order q 2 and not q 3 : 

[P qi , Pq 2 , Pq 3 ] i{QlQ2 + <?3<7l + Q2l3) F ^ Pqi+q2+q 3 

(32) 

This term is reminiscent of the 2-commutator algebra 
Eq.®, and accounts for a possible bidimensional topo- 
logical structure in the 3D insulator. This would be the 
case for a weak 3D Chern insulator, obtained by stacking 
layers of the 2D Chern insulator. This structure remains 
true in all odd dimensions, where the ^-commutator con- 
tains an anisotropic 0(q D_1 ) term in contrast with the 
isotropic (^(q 15 ) term appearing in Ea. (l2Ul) for even di- 
mensions. 

In order to illustrate explicitly the kind of problems 
that arise in odd dimensions, we focus in the following on 
topological insulators in three dimensions. If the Chern- 
Simons density ([29)) is to appear at all in the triple com- 
mutator, this has to be as a 0(q 3 ) term. This term is 
sensitive to the 0(q 2 ) correction to the projected density 
operator, and therefore it matters what kind of regular- 
ization we are considering. In order to be able to make a 
generic statement, we consider in the following the most 
general regularization up to level 0(q 2 ): 



k.n,m 



/lq - E ( 1 - ^ - &<fB MV ) 7 £ f io> <o| 7 r +q 

nm 



(33) 



For instance the parallel transport p q from Eq. (1171) would 
correspond to B^ u — d^Ay — iA^A u . Note that only the 
symmetric part of B^ v is relevant. We computed the 
subdominant term to the 3-commutator (j32[) . The result 
is the sum of two 0(q 3 ) terms, corresponding to the two 
kinds of 3-tensors which are both 0(q 3 ) and antisymmet- 
ric with respect to the three momenta q^. 



r(qi A q 2 A q 3 



E 

k.n.m 



(^AA)„ m7 £t |0) (0| 7 r +qi 



(34) 



and 



k,n,m 



-2^^) nm7 ^|0)(0| 7 ^ qi 



qi+q2+q 3 



(35) 



While the first one e^a^q^q^q^ = (qi Aq 2 Aq 3 )e^ <T 
is familiar from the even dimensional case, the second 
one £ ai a 2 a 3 Qa 1 < la 1 < la2 ^ s on ^y possible in odd dimensions. 
Note that this last tensor being symmetric under ex- 
change of \x and v, it cannot yield the antisymmetric 
Chern-Simons term. At first glance the first term looks 
promising as it contain part of the Chern-Simons form 



6 



(F A A) although the term A A A A A is still missing. 
Before taking a trace though, we need to multiply the 
3-commutator by /9- qi - q2 - q3 , an d one could hope that 
AaAaA will appear. This is not so, as the multiplication 
by p_ qi _ q2 _ q3 yield the additional term 



E 

k.n.m 



(i^4x)wY£ + |0}<0| 7 r (36) 



which kills the term (qi A q 2 A c^)e^ U17 F^ u A a . The final 
result is 

[Plii Pt3-: Pn^P— qi— q 2 — q3 = 

-i E + <#<?3 + ^iX^nrf |0> (0| 7k m 

k,n,m 

+e<. 1 « J «,<^ 1 f^ E (^-^7^10) (0| 7 iT 

k,n,m 

(37) 

where the 3-tensor has the fatal inconvenience to 

be symmetric under ji •<-> v. 



— iDaBpv — id^d v A a (A^d^. + A v d^)A a 
4- F /4 4- F 4 



(38) 



The Chcrn-Simons term being fully antisymmetric, this 
calculation shows explicitly that the Berry curvature can- 
not appear in the algebra of projected density operators 
in three dimensions, no matter what is chosen for the 
regularization B^ v . For instance the parallel transport 
operator p q yields 

T r ([pqi i Pq2 1 Pqs]P-qi-q2-q3 ) = 
+ / rf 2 fcTr(F M „(k)) 



+*w 2 « 3 <« 2 ^ / d2fcTr (^-( k )) + °(q 4 ) 

(39) 

A way to free ourself from this no-go theorem is to 
involve non gauge invariant operators, such as the pure 
translation: 

T q |n> = |n > k-q) (40) 
This can be used to generate the Chern-Simons form as 
1 



P 3 = 



reyjtTr [(p qi p qj (p qfc - T q J- 



g lPq« ^HlKPlj ^q;)(Pqfc ^qjJJP-qi— q 3 — q*, 



(41) 



but the physical picture behind this relation is still un- 
clear. 



CONCLUSION 



We have presented a generalization of the GMP alge- 
bra to D-dimensional topological insulators by general- 
izing the commutator, algebra and Berry phase to their 
higher-dimensional counterparts. This hints at a differ- 
ent group structure from the usual gauge theories. Par- 
ticularly interesting candidates are higher gau ge t heories 
[2ll [27j . In light of this, the recent proposal [241 to de- 
scribe topological insulators by a BF theory [28[ looks 
very promising. In this topological field theory coexist 
two gauge fields. In addition to the usual 1-form A, there 
is an extra D — 1 form B, to which D — 2 dimensional 
membranes can couple. At this level, the even and odd- 
dimensions are fundamentally different - in even dimen- 
sions, the structure factors of the algebra are proportional 
to the D/2'th Chern number, while in odd dimensions 
they are not proportional to the expected Chern-Simons 
form. In two dimensions, the classical limit of the GMP 
algebra is isomorphic to the algebra of area preserving 
diffcomorphisms, and is related to incompressibility. A 
D-algebra on the other hand is related to volume pre- 
serving differomorphisms (29j . Indeed it is a quantiza- 
tion of the classical Nambu-Poisson bracket [23j, which 
is known to be invariant under volume preserving diffco- 
morphisms. It would be interesting to make this con- 
nection more explicit and to understand its link to the 
incompressibility of TIs in higher dimensions. This how- 
ever suggests that the correct "effective" description of 
the higher-dimensional topological insulators is in terms 
of parallel transport not of electrons but of extended ob- 
jects, such as strings in 3 dimensions. We speculate the 
the Chern-simons term could appear when such algebras 
are constructed. 

Note: One day prior to the posting of our paper, a 
related paper[3(| by T. Neupert, L. Santos, S. Ryu, C. 
Chamon, and C. Mudry, containing similar results ap- 
peared on the arxiv. While most of our results are sim- 
ilar, our paper differs in a fundamental aspect. In the 
case of topolo gica l insulators in odd space dimensions, 
the authors of [30j claim that the 9 term appears in the 
odd-dimensional algebras. The basis for this claim, their 
Eq. 2.71, is identical to our equation (jM)) and (133]) when 
we particularize our non-Abelian result to their Abelian 
one. This equation contains the 9 term, although only 
the F A A is present while the A A A A A is missing in 
the non-Abelian Eq. ([53)1. However, Eq. 2.71 of 30] is 
not an algebra, as the right hand side (RHS) does not 
yet contain the p operator and is hence not closed, even 
when the conditions of constant Berry curvature imposed 
in larXivri 202.5188yl are satisfied. In order to close the 
algebra one must recast the RHS of Eq 2.71 of in 
terms of the density, as in our Eq. (I37[) . the 9 term can- 
cels, and is not contain in any of the odd £>-algebras. 
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